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We calculate, in the continuum limit of quenched lattice QCD, the form factor that enters in the 
decay rate of the semileptonic decay B —> Dlv. Making use of the step scaling method (SSM), 
previously introduced to handle two scale problems in lattice QCD, and of flavour twisted boundary 
conditions we extract at finite momentum transfer and at the physical values of the 

heavy quark masses. Our results can be used in order to extract the CKM matrix element Vcb by 
the experimental decay rate without model dependent extrapolations. 



I. INTRODUCTION 



II. FORM FACTORS AND DECAY RATE 



The central goal of flavour physics is the determina- 
tion of the Cabibbo-Kobayashi-Maskawa [3, i quark 
mixing matrix. A set of redundant and precise mea- 
surements can also provide informations about possible 
new physics Q. In turn, precise measurements need an 
adequate theoretical determination of hadron matrix el- 
ements of the weak currents. Non perturbative tools and 
in particular lattice QCD can eventually provide the re- 
quired precision. In this letter we address the determina- 
tion at finite momentum transfer of the form factor that 
enters in the semileptonic heavy meson decay rate. In 
the infinite quark mass limit it is parametrized by the 
Isgur-Wise function [^. In particular, we calculate the 
matrix elements entering the determination of Vcb- The 
calculation that we present misses the effects of dynam- 
ical fermions and is not the final one. Nevertheless, it 
accomplishes for the first time the calculation at finite 
momentum transfer and at the physical values of the 
heavy quark masses, allowing to compare experimental 
data without additional model dependent extrapolations. 

We use a method @ already applied successfully to 
the determination of heavy quark masses and decay con- 
stants [ESUl, called the step scaling method (SSM), and 
special boundary conditions, called flavour twisted 
that shift by an arbitrary amount the discretized set of 
lattice momenta (see also [ 1 ^ [HI). The step scaling 
method allows to reconcile large quark masses with ad- 
equate lattice resolution and large physical volumes; the 
flavour twisted boundary conditions allow to perform a 
calculation at non zero momentum transfer with good ac- 
curacy. We present the main results of our computation 
and stress their phenomenological implications; a more 
detailed discussion on technical aspects will be presented 
elsewhere [l^ . 



The semileptonic decay of a pseudoscalar meson into 
another pseudoscalar meson is mediated by the vector 
part of the weak V — A current. The corresponding ma- 
trix element can be parametrized in terms of two form 
factors. Among possible parametrizations we chose the 
following one 

= {vi -b Vf)f^ h+ + {v, - VfY h_ 




where Mij and Vij = Pij/Mij are the meson masses 
and 4-velocities. The form factors depend upon the 
masses of the initial and final particles and upon w = 

Vf ■ Vi 

h± = = h±{w, Mi, Mf) 

l<w<{Mf + M])/2M,Mf 

In the case where Mi is the B meson mass and Mf is the 
D meson mass the maximum value of w is around 1.6. 
In the infinite mass limit the form factors become 

=^{w) 

M f , Mi 00 

= 0 

where ^(w) is the universal Isgur-Wise function; the con- 
servation of the vector current implies ^(1) = 1. In what 
follows we will And deviations from such a limit that 
will allow a precise discussion on the onset of the HQET 
regime [l^ . 

The differential decay rate for the process B Dlv, 
in the case of massless leptons, is given by 

dw 

\Vcb\^^{MB + MdYMUw^ - 1 ) 3/2 
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The phase space factor — 1)^/^ in the decay rate 
makes its experimental determination harder as w ap- 
proaches 1 and the value at zero recoil is obtained from 
an extrapolation. The function (w), 

= ~ h^f{w) 

^ ^ ^ ' Mf + Mi ^ ' 

is needed in order to extract Vcb by the measurement of 
the decay rate. Previous lattice calculations by the Fer- 
milab collaboration [H,[l3 quote the value of G^^^{w) 
only at zero recoil where it can be extracted with good 
statistical accuracy by using the so called ” double ratio” 
technique. In the following we calculate G^^^{w) in 
the range 1 < w < 1.2 that includes values of w where 
experimental data are available. 



III. LATTICE OBSERVABLES 
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TABLE I: Table of lattice simulations. 



flavour twisted b.c. for the heavy flavours; in particular 
we have used 



We have carried out the calculation within the 0{a) 
improved Schrodinger Functional formalism [I10 with 
T = 2L and vanishing background fields. In order to fix 
the notations, we introduce the following lattice opera- 
tors 

^rs = ^'^CriyhbCi'^) 
y z 



A°^{x) = V'r(a;)757°V's(a;), Prs{x) = llJr{x)j5lps{x) 



leading to 




27rfci 



kiGN 



and ordinary periodic b.c. in the other spatial directions 
and for the light quarks. We work in the Lorentz frame 
in which the parent particle is at rest (pi = 0). In this 
frame w is obtained from the ratio between the energy 
and the mass of the final particle w = Ef/Mf. 

By assuming ground state dominance and by relying 
on the conservation of the vector current, the matrix el- 
ements of can be extracted by considering the ratio 



= {Mf\ \M.)di = 



K^s{x) = tlJr{xh^1ps{x), = tpr{x)^^-f’'lps{x) 



-4°^ (a;) = ^°s(a;) -I- qca 



do 



_^P 

-E r. 



(x) 



2./m;e} 



.Ff_/T/2;0,pf) 

^PU,(TI2- 0, Q>)T)_^{TI2- pf, pf ) 



( 1 ) 



where r and s are flavour indexes while C uud C,' are 
boundary fields at xq = 0 and xq = T respectively. The 
improvement coefficients ca and cy have been taken from 
refs. [iSEIii- We have calculated the following corre- 
lation functions 

^f_./(a;o;Pi,Pf) = '^^{Ou V^f{x) O'ji) 

X 

//(a:o, Pf) = - ^ {OifA]i{x)) 

X 

where i and / refer to the heavy flavours while I to the 
light one. The external momenta have been set by using 



An alternative definition of the matrix elements (D2), 
which reduces to the previous one (Dl) in the limits of 
infinite volume and zero lattice spacing, can be obtained 
by considering 

{v>^y-/ = {Mf\ \Myo2 = 

^ VM-Eff{{T/2,0) .Ff_/T/2;0,pf) 

y/Mjf{{T /2,pf) ^T°^yT/2- 0, 0)T^j^f{T/2- 0, 0) 

( 2 ) 

In eqs. © and © the renormalization factors Zy and 
Za cancel in the ratios together with factors containing 
the improvement coefficients by and 6^. 



3 






FIG. 1: Step scaling functions Lq, I/i), lower plot, 

and Li, L 2 ), upper plot, as functions of the inverse of 

the RGI heavy quark mass rrii of the parent meson measured 
in GeV. The black vertical lines represent the physical values 
of me and mb- The data correspond to the definition Dl. 



By introducing the following ratio 

_ .^;_^(T/2;0,pf) 

(M/I jMf) _ Vw^ - 1 
(M/I yo |M/) “ w + 1 

we can define w, as well as Ef and Mf, entirely in terms 
of three point correlation functions. This definition of 
w is noisier than the one that can be obtained in terms 
of ratios of two point correlation functions; however it 
leads to exact vector current conservation when Mf = Mi 
and reduces the final statistical error on the form factors. 
The two definitions of the matrix elements lead to two 
definitions of 
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FIG. 2: Light quark mass dependence of Lo) iu the 

continuum limit. The different sets of points correspond to 
different values of mi ranging from about ms to about rnsj'^. 
The data correspond to the definition Dl. 



the corresponding RGI quark masses that in the infinite 
volume limit lead to the physical meson spectrum Q ; the 
RGI quark masses are measured by the lattice version of 
the PCAC relation and are not affected by finite volume 
effects (see [H for further details). 

First we compute the observable G^^^{w,Lq) on a 
small volume, Lq, which is chosen to accommodate the 
dynamics of the 5-quark. As in our previous work we 
fixed Lq = 0.4 fm. Then we evaluate a first effect of finite 
volume by evolving the volume from Lq to Li = 0.8 fm 
and computing the ratio 



cr*^^(ui;Lo,Li) 



G^^^{w;Lq) 

G^^D{w,Lq) 



The crucial point is that the step scaling functions are 
calculated by simulating heavy quark masses rrii smaller 
than the 5-quark mass. The physical value {w\ L) is 

obtained by a smooth extrapolation in 1 /m^ that relies on 
the HQET expectations and upon the general idea that 
finite volume effects, measured by the ct’s, are almost 
insensitive to the high energy scale. The final result is 
obtained by further evolving the volume from L\ to L 2 = 
1.2 fm, according to 






r 



1 '^ wr-1 

M,(l + r) ^ {Vy^f J 



Ml 

M, 



( 3 ) 



The last equation is not defined at w = 1; this is due 
to the second term in parenthesis that we extrapolate at 
zero recoil before calculating G'^^^{w = 1). 



G^^°{w-L2) = 

= G^^°{w-Lq) a^^°{w-LQ,LQ) 

Physical values require also usual continuum and chiral 
extrapolations. 



V. RESULTS 



IV. THE STEP SCALING METHOD 

The SSM has been introduced to cope with two-scale 
problems in lattice QCD. In the calculation of heavy-light 
meson properties the two scales are the mass of the heavy 
quarks (5,c) and the mass of the light quarks (u,d,s). 
Here we consider the form factor G*^-^ as a function of 
rc, the volume, L^, and identify heavy meson states by 



In figure [T] we can test the validity of the SSM. The 
step scaling functions are almost insensitive to the heavy 
quark mass rrii of the parent meson for values larger 
than rric. Moreover, finite volume effects are already very 
small at L = 0.8 fm, in particular at zero recoil. 

We present results already extrapolated to the chiral 
and continuum limits. We have simulated three different 
lattices for Lq and two for the other volumes (see table |T| 
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FIG. 3: Comparison of the two definitions of L) at 

I/O = 0.4 fm (upper plot) and at L 2 = 1.2 fm (lower plot). 
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TABLE II: Final result. Form factor in the con- 

tinuum and infinite volume limits. As a comparison we quote 
also the results of previous lattice calculations by the Fermilab 
collaboration. 

observing small 0{a) effects on the form factor and on 
the step scaling functions. Further details on the contin- 
uum extrapolations will be given in ref. [l^ . Concerning 
the chiral behaviour, our results do not show any sizable 
dependence upon the light quark mass, as shown in fig- 
ure [2] for The same feature is observed 

for all the combinations of heavy quark masses and on 
all the volumes that we have simulated. Nevertheless we 
make a linear extrapolation to reach the chiral limit; the 
resulting error largely accounts for the systematics due 
to these extrapolations. 

As discussed in sec. Ell we used two definitions of the 
form factor that, at finite volume, differ by the finite vol- 
ume effects. A check of the convergence of our SSM can 
be obtained by comparing these two definitions on the 
smallest volume (figure |3] upper plot) and on the largest 
one (figure [3] lower plot). We see that, while the small 
volume results differ, the final ones converge to a com- 
mon value giving us confidence of a correct accounting of 
finite volume effects. 




FIG. 4: Comparison of \Vcb\ G^~'^{w) with available exper- 
imental data. The plot has been done by normalizing our 
lattice data with the value of Vcb extracted at ui = 1.2. 



VI. CONCLUDING REMARKS 

In table El we quote our final results obtained by av- 
eraging over the two definitions and by combining in 
quadrature statistical errors with the systematic ones due 
to the small residual dispersion between D1 and D2. As 
a comparison we show previous lattice results obtained 
by the Fermilab collaboration at zero recoil. 

The existence of predictions up to w ~ 1.2 and physical 
b and c quark masses allows a direct comparison with 
experimental data, as shown in figure E The comparison 
has been done by extracting the value of Vet by the ratio 
of the experimental and lattice data at w = 1.2; as an 
indication, we get Kh = 3.84(9) (42) x 10“^, where the 
first error is from our lattice result, = 1.2) = 

0.853(21), and the second from the experimental decay 
rate, \Vcb\G^^^{w = 1.2) = 0.0327(35), as deduced from 
the plots of refs. [10, [2l|, . 

The extension of this calculation to the unquenched 
case does not present problems of principle. The recur- 
sive matching process can be extended to the sea quark 
masses that, alternatively, can be kept to their physical 
values if the Schrodinger Functional formalism is used. 
Moreover, flavour twisted boundary conditions can be 
used for heavy valence quarks also in the Nf = S un- 
quenched theory. The real case will further differ by the 
heavy flavour determinants that can be computed per- 
turbatively. 
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